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SUPERSONIC FLUTTER OF A THERMALLY STRESSED FLAT PANEL 

WITH UNIFORM EDGE LOADS* 

By Harry G. Schaeffer and Walter L. Heard, Jr. 
Langley Research Center 


SUMMARY 

The effect of a nonuniform, self- equilibrating stress state associated 
with a parabolic temperature distribution is included in the flutter analysis 
of a rectangular, simply supported panel subjected to airflow over one surface 
and uniform edge loads. Linearized, static, two-dimensional aerodynamics is 
used in conjunction with thin-plate theory to formulate the problem. Numerical 
results for a square panel show that nonuniform self -equilibrating stresses are 
as important as uniform stresses in affecting flutter behavior. 


INTRODUCTION 



Dynamic instability of thin platelike structures in a gas flow is a prob- 
lem which has received considerable attention in recent years. As a result of 
the high velocities attainable by modern aircraft, severe aerodynamic loading 
of exterior skin panels frequently occurs and can lead to panel flutter. (See, 
for example, ref. 1.) In the treatment of panel-flutter problems, theoretical 
investigators have taken into account such effects as panel configuration, flow 
orientation, boundary conditions, and uniform midplane stress. However, effects 
of thermally induced, self-equilibrating, midplane stress have not been 
examined. 

It has already been shown in reference 2 that a system of uniform stresses 
can greatly reduce the flutter speed of an unbuckled panel. The purpose of the 
present paper is to investigate the effects of a nonuniform, self-equilibrating, 
stress state on the flutter behavior of a flat, uniformly stressed panel. Since 
a nonlinear temperature distribution could easily result from aerodynamic 
effects, it is assumed for the present analysis that the nonuniform stress state 
is thermally induced by a parabolic temperature distribution which is constant 
through the thickness. However, the results are valid for any similar system 
of self -equilibrating midplane stresses. 


*Part of the information presented herein was offered as a thesis by the 
second author in partial fulfillment of the requirements for the degree of 
Master of Science in Engineering Mechanics, Virginia Polytechnic Institute, 
Blacksburg, Virgina, December 1964. 


For the purpose of this analysis the panel is considered to be simply 
supported against lateral deflections, but unrestrained to inplane displace- 
ments at the boundaries. In order to simplify the analysis the aerodynamic 
forces are assumed to be those resulting from Ackeret theory. The analysis 
proceeds by using the Galerkin technique in conjunction with an assumed stress 
function to determine an approximate midplane stress distribution corresponding 
to the assumed temperature variation and boundary conditions. This thermal 
stress distribution is then substituted into the linear partial differential 
equation governing lateral deflections of thin plates in combined bending and 
compression. The Galerkin procedure is used again to obtain approximate solu- 
tions for the flutter behavior, and numerical results are presented for a 
square panel subjected to various combinations of uniform loading applied at 
the boundaries-. 


SYMBOLS 

amn amplitude coefficients (see eq. (ll)) 

a panel length 

b panel width 



E Young ' s modulus 
h panel thickness 
i,j,m,n,r, s integers 


k 

o 7ha 2 
k^ = Cu 

jt^D 

(I 1 )rs^ ( I 2) rs ^ ( I 5)rs integrals defined by equations ( 13 b) to (ljjd) 
K r s = Mrs + p rs “ ^Qrs “ 
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aerodynamic loading coefficient 


j rs 


l aerodynamic pressure load per unit area 

M Mach number 

M rs coefficient defined by equation (A5) 

N x , Ny , Nxy midplane stress resultants 

Nx 0 ,% 0 midplane uniform stress resultants 

Nx^Ny^N-jjy^ midplane thermal stress resultants 


r rs 

Qrs 

R 

Rv_ = 


coefficient defined by equation (A 6) 
dynamic pressure 

coefficient defined by equation (A7) 
total number of terms in flow direction 


N x a 2 
it 2 D 


%o 

S 

t 

T 


= %o t 

3t 2 D 


total number of terms in cross-flow direction 
time 

temperature 

maximum amplitude of temperature distribution 


w 

x,y,z 

a 

7 


fl 



lateral deflection of panel 
Cartesian coordinates of panel 
coefficient of thermal expansion 
mass density of panel 
when i =5 1 
when i / j 





A C r 

9 


dynamic-pressure parameter, 

critical value of A 
Poisson's ratio 
Airy stress function 

thermal stress parameter, 


2qa- 


D\Zm2~ 


aEha^ATi 
it 2 . D 



ft generalized eigenvalue 

cd frequency 

Notation 



^ - ( ),xxxx + 2 ( ),xxyy + ( 

When subscripts follow a 
of the principal symbol with respect to the subscripts. 


),yyyy 

comma, they indicate partial differentiation 


STATEMENT OF PROBLEM AND ANALYSIS 


A flat, rectangular panel of length a, width b, and uniform thickness h, 
as shown in figure 1, is considered in this analysis. The panel is simply sup- 
ported on all edges, exposed to supersonic airflow at Mach number M over one 
surface, and subjected to a nonlinear temperature distribution T = T(x,y) in 
the middle plane. The panel edges are considered to be unrestrained to thermal 
expansion; however, uniform loading, Nx 0 and Ny Q , which is independent of 

the temperature distribution may be applied at the boundary. 


Basic Equations 

In this analysis thin-plate theory including the effects of inplane loads 
is used. The basic equilibrium equations for this case are 


iaAw + 


N x w ,xx 


2N xy w ,xy + 


N y w ,yy 


+ 7hw 


( 1 ) 
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N x,x + N xy,y 
N y,y + N xy,x 



( 2 ) 


It is assumed that N x and Ny are positive in compression and that all mate 
rial properties are to be evaluated at the average panel temperature. 


An approximate expression for the aerodynamic pressure load Z is obtained 
by use of linearized, static, two-dimensional aerodynamics (Ackeret theory) . 

This approximation not only greatly simplifies the analysis, but as shown in 
reference 3, yields essentially the same flutter results as exact aerodynamics 
at M ^ 1.3 for uniformly stressed panels, except when the stress level causes 
two modes of oscillation to have the same, or nearly the same, frequency. Thus, 
if care is taken to avoid these stress levels, Ackeret theory should also be 
applicable to the nonuniform stress problem of the present analysis, so that 


Z 


-2q 

\/m2 - 1 


(3) 


The midplane force intensity terms in equations (l) and (2) are written as 
the stun of those induced by applied, uniform, normal forces at the boundary 
(subscript o) and those resulting from the nonlinear temperature distribution 
(subscript T) as follows: 


N x = N Xq + 


N y = N y 0 + %t 




(M 


N xy ~ N xy T 


J 


The thermally induced stresses are determined in terms of a stress func- 
tion <p = cp(x,y) defined by the following relations: 


Nx t = y,yy ^ 

N y T = 9, xx > 
N xy T = -<P,X£ 


(5) 


This stress function identically satisfies the equilibrium conditions given by 
equations (2). For compatibility of inplane strains, the stress function must 
satisfy the following partial differential equation (ref. 4): 

V^cp = ccEhV 3 ! 1 (6) 


where a is the coefficient of thermal expansion at the average plate tempera- 
ture and the temperature T is a function of x and y. Note that the sign 
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of the right-hand side is opposite to that normally encountered, because direct 
stresses are positive in compression. 

The panel boundary conditions have been chosen in such a way that they can 
be satisfied by using rather simple functions. Thus, the resulting flutter 
analysis is greatly simplified, but still predicts trends which one may expect 
to be valid for the actual panel. 


The condition that the panel be free from thermally induced normal and 
shear stresses on the boundaries requires that the stress function satisfy the 
following boundary conditions: 

<p(0,y) = <p(a,y) = cp(x,0) = q>(x,b) = 0 

<P,x(°>y) = <P,xKy) = <P,y(x,0) = cp,y(x,b) = 



The condition that the panel boundaries be simply supported requires that 
the displacement function satisfy the following boundary conditions: 


v (0,y, t ) = w(a,y,t) = w(x, 0,t ) = w(x,b,t) = 0 

(0,y,t ) = w^ xx (a,y,t ) = w^yy(x,0,t) = w^yy(x,b,t) = 0 1 


( 8 ) 


w 


XX' 


The analysis which follows consists of two parts. First, the coefficient of an 
assumed stress function cp is determined approximately by using the Galerkin 
technique; then, the stability determinant is determined from the deflection 
equation, again by using the Galerkin technique. 


Solution of Stress Function Equation 

The temperature distribution, considered in this problem, varies paraboli- 
cally in the flow and cross-flow directions and may be represented by 

T - 1^(|)(X - i)(£)(l - f) 

where ATp is the maximum amplitude of the temperature distribution. 

In order to simplify the subsequent flutter analysis, it is assumed that 
an approximate stress function is of the following form: 

9 - 0«Eha 2 ^ 1 (|) 2 (| - l) 2 (j0 2 (| - 1 f (9) 

This approximation of the stress function satisfies all the boundary conditions 
given by equations ( 7 ), but does not satisfy equation (6). The best fit is made 
by applying the Galerkin technique to equation (6) to evaluate the constant C. 
The result is 


6 


c 



( 10 ) 


The accuracy of equation ( 9 ) in approximating the true stress state associated 
with the temperature distribution will be assessed subsequently. 


Solution of Differential Equation 

It is assumed that the solution of equation (l) can be represented as 
follows: 


00 00 

w(x,y,t) = Re ^ ^ ^n sin TT sin IT 

m=l n=l 


where the frequency o> is, in general, complex. Each term of the assumed 
solution satisfies the given boundary conditions on w. 

After substituting equations ( 3 ), (^), (5)* (10), and (ll) into equa- 
tion (l) and applying the Galerkin procedure, the following set of equations is 
obtained for the amplitude coefficients a ^ n : 


a\ 2 o2 _ v 2l 


r 2 + ^ s 2j . RxQr 2 _ Ryo ^ s 2 _ k L + * Lrs 


Jt_ £ 

*2d b 


( X l) 


rs 


+ 2 tfe): + 


/rs 


(b) ( I 3)rsJ 


( 12 ) 


where 


Lrs * L L ll- cos ^ sin ^ sin t sin t ^ 


(i» 


m=l n=l 


( I l)rs - flfl l X sin =2 sin ZS- sin HZ sin SS to dy (15b) 

w U -i -i 


m=_L a=x 


f -r . N 

V^rs 


r h r & v v 


~ J J l-j L mnamnC( W cos TT sin TT cos IT sin IT ** dy 


o u o 


m=l n=l 
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00 00 

- J b L I I n 2 *WP,xx sin 5KS sin 22£ sin 22 sin £2 dx dy ( 13 d) 
0 0 m=l n=l 

r = 2^ 3 ^ • • • 

S 5 1; 2 ; 3, • • • 


and 


2 7ha\)2 

k = ls~ 


R 


Nx 0 a 2 

Xo = ~7* 


R 


N yo a2 

‘ y ° = « 2 d 


The integrals (dl) rs > (^2) rs ^ and (•'■3) r s arise due to the nonuniform 

stress distribution and the integral L rs which arises due to aerodynamic 

loading are evaluated in appendix A. By using the results of appendix A, equa- 
tion (12) becomes 




- R x 0 r2 " 



^rs + ^(tf* K rs = 0 


(14) 


where 


r — ly 2, 3> • • • R 
s = 1, 3, 5 , . . . 2S-1 

and 


* = 


aEha 2 ^-^ 

3t%) 
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The coefficients L rs , M rs , P rs , an< ^ Qrs are defined by equations (Al), 

(A 5 ), (a 6 ), and (A7), respectively. 

The parameter \|r is analogous to the nondimens ional, midplane, uniform 
stress parameters Rx 0 and Ry Q and may be thought of as a thermal stress 

parameter. The range of the indices r and s is R and 2S-1, respectively, 
where R is the total number of sine terms in the flow direction and S is 
the total number of sine terms in the cross-flow direction. Only the sine 
terms corresponding to an odd number of half waves in the cross-flow direction 
need to be considered, since, for the symmetric temperature distribution used, 
even and odd cross-flow terms uncouple and the odd terms lead to the lowest 
(critical) flutter speed. Thus, equations (14) represent RS (R-times-S) 
linear, homogeneous, algebraic equations for the RS unknown amplitude 
coefficients a rs . In order that nontrivial solutions of the system of equa- 
tions (lk) exist, it is necessary that the determinant of coefficients vanish. 
Thus, the stability criterion may be written in the following form: 

det(Aij - =0 ( 15 ) 

where ft is the eigenvalue. 

Since the problem is one of determining the stability of a given form of 
solution, it is most natural and advantageous to associate the eigenvalue with 
the frequency parameter k^. Then, for the nonconservative problem, the system 
is dynamically unstable when the eigenvalue ft becomes complex. This implies, 
according to equation (ll), that the system diverges in an oscillating fashion. 
The problem now is to determine the relationship that must exist between the 
parameters a/b, R Xo > Ry 0 , k^, and ^ or( ^ er that the system be stable 

in the sense described. 

For the case where R = 2 and S = 1, that is, 
w = ^a]_p sin + a 2 q sin ^~sjsin 

equation ( 15 ) simplifies to the point where it is possible to determine an 
analytical expression for values of A corresponding to harmonic motion of the 
panel. This result is 


9 



(16) 




0.4c 








The critical value of A can be found by maximizing A, expression (l6), 
with respect to k2. Solving for k^ and substituting the result into equa- 
tion (l6) gives 


Xcr 




+ 


O.lgC 

*2 



(17) 


where A cr is the critical value of the dynamic-pressure parameter which leads 
to flutter. 


For higher ordered modal solutions, the algebraic difficulties in obtaining 
expressions similar to equations (l 6 ) and ( 17 ) are formidable, and recourse is 
made to an iterative procedure using a digital computer. For the computer solu- 
tion the eigenvalues are calculated, and the lowest value of A for which two 

of the eigenvalues coalesce is sought. The approximate modal solutions con- 
sidered are for the cases where R = 4, S = 1; R = 6 , S = lj R = 6 , S = 2; 

and R = 6 , S = 3- The coefficients associated with the latter case have been 

calculated from equations (l4), and the general matrix equation is presented in 
appendix B. The coefficients associated with the other cases can be determined 
by deleting appropriate rows and columns. 


DISCUSSION OF RESULTS 


Thermal Stress Distribution 

An approximate stress distribution has been assumed in order to simplify 
the subsequent dynamic analysis. Since the purpose of this paper is to deter- 
mine the effect of the stress state associated with a nonlinear temperature dis- 
tribution on the flutter characteristics, it is important to consider the accu- 
racy with which the stresses obtained from the assumed stress function approxi- 
mate the true stress state. The true stress state associated with the parabolic 
temperature distribution has been determined by the finite-difference method of 
reference 5 * 

Comparisons between the one-term approximation and a finite-difference 
solution are presented in figures 2 and 3 . Figure 2 shows the variation of 
Nxrp with x for specified values of y, and figure 3 shows the variation of 

Nxyrp with x for specified values of y. The normal stress determined by the 
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one-term approximation and that determined by the finite-difference solution 
agree quite well (fig. 2). In the case of the shear stress distribution 
(fig. 3 ), the agreement in magnitude is not as close as that for the direct 
stresses; however, it appears that the agreement is sufficiently close so that 
the overall results can be interpreted with confidence. 


Effect of Temperature Distribution on Flutter Boundary 

The effect of the parabolic temperature distribution on the flutter 
behavior of a square panel free of uniform midplane stress is shown in figure k. 
The results are presented for a 6-by-3 (R = 6 ; S = 3) term solution. The inter- 
section of the flutter boundary with the A-axis is the value of A cr associated 

with the unstressed, unheated panel. This value decreases as ijr increases 
(ATp increases) until the flutter boundary becomes tangent to the curve labeled 

"thermal buckling loop." 

The three regions shown in figure 4 are characterized by the value of the 
frequency parameter squared (k^) . In the region labeled "panel flat, no flut- 
ter, " k^ is real and positive; hence, cn is real. In the region labeled 
"flutter, " k^ is complex; thus, u) is complex and at least one root will lead 
to oscillating, divergent panel motion. In the region labeled "panel buckled, 
no flutter," k^ is negative; thus, to is pure imaginary and the panel is 
statically unstable. 

These three regions are separated by two boundaries. The first is the 
buckling loop which is the locus of points for which k^ = 0. The second is 
the flutter boundary which is the locus of points at which two frequencies 
coalesce. The flutter boundary is terminated at its point of tangency with the 
static buckling loop, since this point represents the limit of linear plate 
theory. This transition point is of considerable importance, because experience 
has shown that it represents the lowest value of A cr associated with this 

panel configuration. 

Figure 4 shows that the stress distribution associated with a parabolic 
temperature distribution can cause a 6 l-percent reduction in the value of Acr 
associated with an unheated, unstressed panel. In order to give an indication 
of the magnitude of ATq 'which causes an effect of this order, ATq is com- 
puted for a representative square aluminum panel with ^ = 300. The results 

show that a temperature difference of only 27° F between the center and the 
edges of the panel causes the 6 l-percent reduction in Acr- 

The effect of the parabolic temperature distribution on the flutter bound- 
ary associated with uniform compressive loads in the flow direction is shown in 
figure 5 for the 6 - by-3 term analysis. The curve labeled "\|r = 0" is essentially 
the exact solution established in reference 2 and shows how the critical value 
of the dynamic -pres sure parameter varies with applied uniform loading in Lhe 
flow direction. Increasing values of t lower the flutter boundary. 
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The regions of stability and instability are not shown in figure 5 but are 
similar to those shown in figure however, in figure 5 the dynamic-pressure 
parameter is plotted against Rx 0 instead of >| r. For every curve represented 

by some constant value of \|r, there is an associated buckling loop. The dashed 
line represents the locus of points where the flutter boundary becomes tangent 
to its corresponding buckling loop. 


The decrease in A cr due to a parabolic temperature distribution is of the 
same order of magnitude as the decrease in Acr due to a uniform compressive 
load in the x-direction. The comparable reductions in A C r due to the two dif- 
ferent parameters indicate that the effect of the nonuniform temperature dis- 
tribution is as important as that of uniform compressive stress. Since the 
nonuniform stress distribution could have resulted from causes other than a 
temperature distribution, the previous statement can be generalized to say that 
any system of nonuniform self-equilibrating stresses may have a significant 
effect on panel flutter behavior. 


Reference 2 shows that uniform compressive loads applied perpendicular to 
the direction of airflow have virtually no effect on the flutter of a uniformly 
stressed panel. It is of interest, therefore, to investigate the effect of 
Ry 0 for the present case. In figure 6 the relation between Acr an d Rx 0 

presented for various values of the stress ratio R y 0 y /^x 0 * Th e differences in 
the flutter boundaries due to Ry Q are not significant and cannot be plotted 

on the scale presented. However, the transition point between the flat-panel 
and the buckled-panel flutter boundaries is significantly affected as indicated 
by the loci of termination points of the flat-panel flutter boundaries. Thus, 
for the given range of parameters and the assumed nonuniform stress distribu- 
tion, the results of figure 5 ma y be used to predict the flutter behavior of a 
square panel with good accuracy for any combination of Rx 0 > Ry 0 > anc ^ ^ which 


does not cause the panel to buckle. Only the terminal points need to be 
located. The loci of these terminal points for several values of the stress 
ratio are shown in figure 6 . The curves apply for the stress-ratio range 


0 < 


R- 




R- 


*0 


< 1 


only, since large positive values and large negative values of the 


ratio may yield different results because of changes in buckling mode shapes in 

R y 

the cross-flow direction. When — — = 1 , the locus of termination points of 

Ry- _ 

A 0 

the flat-panel flutter boundary is essentially a horizontal line. Thus for 
this stress ratio, the theory predicts that the transition point is essentially 
independent of the relative magnitudes of applied uniform stress and thermal 
stress. 


Since the results presented in this paper were obtained by applying the 
Galerkin technique, there may be some question as to convergence toward the 
exact solution of the problem. In order to investigate the question of con- 
vergence, results for various approximate solutions are presented in figure 7 - 
The flutter boundaries of A cr as a function of R Xq for values of 
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t of 10 , 20 , 30 , and 40 are presented in figure 7 for the 2 -by-l, 4-by-l, 

6 -by-l, 6-by-2, and 6 -by - 3 term analyses. The solid curves represent results 
of the 6 -by -2 and 6 -by -3 term analyses whereas the other curves represent the 
results of the lower order solutions. 

An examination of figure 7 shows that the flutter boundary is significantly 
altered when the number of terms in the flow direction is increased from two to 
four. However, a further increase to six terms exhibits very little effectj 
thus, the solution is assumed to be converged for flow- wise terms. The require- 
ment of relatively few flow-wise terms for convergence is not unexpected, since 
for this range of parameters a four-term solution gave a close approximation to 
the exact solution in reference 2 where is zero. 

Inclusion of two cross-flow terms in the analysis instead of only one 
exhibits only a very slight lowering of the flutter boundary as may be seen by 
comparison of the 6 -by-l and 6-by-2 solutions. (For the scale shown, the 
6 -by-l and 4-by-l results are coincident except in fig. 7 (a) where the 6 -by-l 
solution gives the higher values of A C r*) Inclusion of three cross-flow terms 
gives essentially the same results as the 6 -by -2 solution and indicates that 
the 6 -by -2 solution is converged for the range of <|r considered. 


CONCLUDING REMARKS 


Nonuniform stresses associated with a parabolic temperature distribution 
are shown to be as important in affecting panel flutter behavior as uniform 
loading applied in the direction of airflow. It Is also shown that although 
uniform loading applied perpendicular to the direction of airflow has virtually 
no effect on the flutter boundaries for values of the stress ratio Ry 0 /Rx 0 

from 0 to 1, the transition point between the flat-panel flutter boundary and 
the buckled-panel flutter boundary is significantly affected. As this stress 
ratio approaches the value of 1, the transition point becomes essentially inde- 
pendent of relative magnitudes of applied uniform loading and thermal stress. 

Because of the effects of the nonuniform stresses shown herein, it appears 
that thermally induced nonuniform stresses should be considered in the correla- 
tion of experimental flutter results with theory. In fact, consideration should 
be given to any system of nonuniform, self-equilibrating stresses, even those 
resulting from causes other than a nonlinear temperature distribution. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 12, 1965 . 
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APPENDIX A 


LOAD COEFFICIENTS 


Evaluation of the integral resulting from the aerodynamic loading 
(eq. (13a)) gives the following expression: 


R 



m/r 


r-4 1 - ( - 1)r j 

r^ - m^ 


(Al) 


where the notation m £ r indicates that the terms in the sum associated with 
m = r are to be omitted • 

The integrals (Ii) rs * ( I 2) rs > and ( I 3) r s given by equations (13b), 

(13c), and (13d), respectively, are evaluated by making use of the specific 
definition of the stress function given by equation (9). The result is 


Mrs 

(■^2)rs 

Mrs 


2CaEha : 7>T 1 

it^b 

It-CoEha 2 ^ 


M- 


TS 


Or 


TS 


2CaEhabAT^ 


■ rs 


7t c 


(A2) 

(A3) 

(Alv) 


where 


R 2S- 

I l 

m=l n-1, 3 
m^r n/-s 


1 Mnst-l)^ 8 + l] [(-l) m+r + l] [(. + - (m - r)J 


a mn 


12 

^nt) 2 

tfi 

60 




2S-: 

z 

n=l, 3 
n/s 


:A(n - _ s) 2 (n + s) 2 (m + r^m - r)^ 

1 .n+s 1 

ns[(-l) lja rn 


(n - s) 2 (n + s) 2 



rn/r 



m£[(-l) m+r + l] [(m + r) 4 - (> - r a ffls 
(m + r)^(m - r)^ 


+ 



a rs 


(A5) 
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APPENDIX A 


R 2 S -1 

Z I 

m=l n=l ,3 
mfr n/s 


a mn\ 


jWw[(-l) m+r + l] [(-1)“ +B + l] [(n + s) 4 - (n - S ) 
:A(n - s )^( n + s) i+ (m + r) 2 (m - r) 2 


12 

(sit ) 2 + 

3 

R H N m+r 

Y 1 mr[(-l) +lja ms 


* 2 

60 

4 (sxc ) 2 

(”> - r) 2 (m + r ) 2 






m^r 


9 

2 S -1 

V 

n 2 [l + (-l) n+B ] [(„ + s)* - (n - s)* 

|_ a 

t 2 i^ 

n- 1,3 


(n + s)^(n - s ) 4 

^rn 



n/s 




3 c- 

A 2 

i + 

3 



2 V 


60 


a rs 



(A 6) 


H 2 S -1 


«»-Z Z ? -[<-«“♦ i [<-»■”* i 


EQr=l tt. 1,3 

m^r n^s 


f (s - np + (s + np 

(r - mp + (r + mp 1 

^ (s - np(s + np 

(r + mp(r - mp J 


W^s 2 




m=l 

m^r 

2 S -1 


(r - m)^ + (r + m)^ 


(r - m)^(r + mp 


a lBS 


*Z>>“*3 


n*l ,3 

n/s 

16A 2 s 2 TS 


(s - n)^ + (s + n)^ 


(s - n)^(s + n)-^ 


a rn 


(A7) 


As in equation (Al), the notation m f r, n / s indicates that the terms in 
the sums associated with m = r and n = s are to be omitted. 
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ONE-TERM APPROXIMATION 
FINITE-DIFFERENCE SOLUTION 



Figure 3.- Variation of shear stress resultant N X y T with x. T = ltoTjLU - -Mu - ^ = 1; antisymmetric atxrut 
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Figure 4.- Thermal flutter boundary for square pane! with no externally applied forces. 
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